Improved preconditioned conjugate gradient algorithm and application in 3D inversion of gravity-gradiometry data   总被引：1，自引：0，他引：1  

Tai-Han Wang  Da-Nian Huang  Guo-Qing Ma  Zhao-Hai Meng  Ye Li 《应用地球物理》2017,14(2):301-313

With the continuous development of full tensor gradiometer (FTG) measurement techniques, three-dimensional (3D) inversion of FTG data is becoming increasingly used in oil and gas exploration. In the fast processing and interpretation of large-scale high-precision data, the use of the graphics processing unit process unit (GPU) and preconditioning methods are very important in the data inversion. In this paper, an improved preconditioned conjugate gradient algorithm is proposed by combining the symmetric successive over-relaxation (SSOR) technique and the incomplete Choleksy decomposition conjugate gradient algorithm (ICCG). Since preparing the preconditioner requires extra time, a parallel implement based on GPU is proposed. The improved method is then applied in the inversion of noisecontaminated synthetic data to prove its adaptability in the inversion of 3D FTG data. Results show that the parallel SSOR-ICCG algorithm based on NVIDIA Tesla C2050 GPU achieves a speedup of approximately 25 times that of a serial program using a 2.0 GHz Central Processing Unit (CPU). Real airborne gravity-gradiometry data from Vinton salt dome (southwest Louisiana, USA) are also considered. Good results are obtained, which verifies the efficiency and feasibility of the proposed parallel method in fast inversion of 3D FTG data.  相似文献   

Synchronous pattern of moisture availability on the southern Tibetan Plateau since 17.5 cal. ka BP – the Tangra Yumco lake sediment record     

Marieke Ahlborn  Torsten Haberzettl  Junbo Wang  Karoline Henkel  Thomas Kasper  Gerhard Daut  Liping Zhu  Roland Mäusbacher 《Boreas: An International Journal of Quaternary Research》2017,46(2):229-241

A possible asynchronicity of the spatial and temporal moisture availability on the Tibetan Plateau has been a controversial subject of discussion in recent years. Here we present the first attempt to systematically investigate possible spatial and temporal variations in moisture availability by examining two lakes, Tangra Yumco and Nam Co, on an east–west transect on the southern Tibetan Plateau using identical proxies for palaeoenvironmental reconstruction. In this study, an independent record from Tangra Yumco was analysed applying a multi‐proxy approach to reconstruct variations in moisture availability since the Lateglacial. Results were subsequently compared with previously published records from Nam Co and additional records from Tso Moriri (northwestern Himalaya) and Naleng Co (southeastern Tibetan Plateau). Our results show that Tangra Yumco was at least partially ice covered prior to 17.1 cal. ka BP. A temperature rise after 17.1 cal. ka BP probably resulted in thawing of the permafrost. At 16.0 cal. ka BP moisture availability increased, representing an initial monsoonal intensification. Warmer conditions between 13.0 and 12.4 cal. ka BP and cooler conditions between 12.4 cal. ka BP and the onset of the Holocene reflect the Bølling‐Allerød and Younger Dryas. At the onset of the Holocene moisture availability rapidly increased, with moisture highest prior to 8.5 cal. ka BP when temperatures were also highest. After 8.5 cal. ka BP the moisture availability gradually decreased and showed only minor amplitude variations. These findings are consistent with the records from large lakes like Nam Co, Tso Moriri, and Naleng Co, revealing a synchronous pattern of moisture availability on the southern Tibetan Plateau.  相似文献   

Elastic lateral dynamic impedance functions for a rigid cylindrical shell type foundation          下载免费PDF全文

R. He  R.Y.S. Pak  L.Z. Wang 《国际地质力学数值与分析法杂志》2017,41(4):508-526

Elastic lateral dynamic impedance functions are defined as the ratio of the lateral dynamic force/moment to the corresponding lateral displacement/rotation at the top ending of a foundation at very small strains. Elastic lateral dynamic impedance functions have a defining influence on the natural frequencies of offshore wind turbines supported on cylindrical shell type foundations, such as suction caissons, bucket foundations, and monopiles. This paper considers the coupled horizontal and rocking vibration of a cylindrical shell type foundation embedded in a fully saturated poroelastic seabed in contact with a seawater half‐space. The formulation of the coupled seawater–shell–seabed vibration problem is simplified by treating the shell as a rigid one. The rigid shell vibration problem is approached by the integral equation method using ring‐load Green's functions for a layered seawater‐seabed half‐space. By considering the boundary conditions at the shell–soil interface, the shell vibration problem is reduced to Fredholm integral equations. Through an analysis of the corresponding Cauchy singular equations, the intrinsic singular characteristics of the problem are rendered explicit. With the singularities incorporated into the solution representation, an effective numerical method involving Gauss–Chebyshev method is developed for the governing Fredholm equations. Selected numerical results for the dynamic contact load distributions, displacements of the shell, and lateral dynamic impedance functions are examined for different shell length–radius ratio, poroelastic materials, and frequencies of excitation. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Potential wind erosion rate response to climate and land‐use changes in the watershed of the Ningxia–Inner Mongolia reach of the Yellow River,China, 1986–2013          下载免费PDF全文

Heqiang Du  Tao Wang  Xian Xue 《地球表面变化过程与地形》2017,42(13):1923-1937

Climate and land‐use changes could strongly affect wind erosion and in turn cause a series of environmental problems. Thus, the objective of this study was to assess potential wind erosion rate (PWER) response to climate and land‐use changes in the watershed of the Ningxia–Inner Mongolia Reach of the Yellow River (NIMRYR), China. The watershed of NIMRYR suffers from serious wind erosion hazards, and over recent decades, wind erosion intensity and distribution has changed, following climate and land‐use changes. To understand these processes in the NIMRYR watershed, the Integrated Wind Erosion Modelling System (IWEMS) and the Revised Wind Erosion Equation (RWEQ) were used to calculate the PWER under different climate conditions and land‐use scenarios, and to assess the influences of climate and land‐use changes on the PWER. The results show the PWER in the whole watershed had a significant declining trend from 1986 to 2013. The results of the relationship among PWER, climate change, and land‐use changes showed that climate change was the dominant control on the PWER change in this watershed. Compared to the period 1986–1995, the average PWER decreased 23.32% and 64.98% as a result of climate change in the periods 1996–2005 and 2006–2013, respectively. In contrast with climate change, the effects of land‐use changes on the average PWER were much lower, and represented a change in PWER of less than 3.3% across the whole watershed. The study method we used could provide some valuable reference for wind erosion modelling, and the research results should help climate and land‐use researchers to develop strategies to reduce wind erosion. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

Comparison of the effects of unidirectional and sign‐alternating temperature gradients on the sintering of ice spheres          下载免费PDF全文

Xuan Wang  Ian Baker 《水文研究》2017,31(4):871-879

To simplify the complex snow structures that occur in nature, polycrystalline ice spheres were produced and arranged vertically to model the sintering process. By controlling the temperatures on both the top and bottom of the ice sphere array, the effect of upward and downward vapor transfers was examined. The evolution of the neck areas between ice spheres was observed using X‐ray computed microtomography. As frequently observed under the basal part of a snow layer and previous experiments of snow temperature gradient metamorphism, depth hoar structures were formed along neck areas and their formation was found to be directly related to the vapor transfer direction. To model the temperature gradient inversion that can be induced in nature by daily cycles of radiative heating and cooling, we also performed sign‐alternating temperature gradient experiments on the ice sphere arrays. The morphological evolution of the neck and the associated vapor transfer were examined through image analysis and 2D modeling. The final microstructures of the neck area turned out to be a symmetrical distribution of ice protrusions bridging neighboring ice spheres.  相似文献   

Three‐dimensional shakedown solutions for anisotropic cohesive‐frictional materials under moving surface loads     

Juan Wang  Hai‐Sui Yu 《国际地质力学数值与分析法杂志》2014,38(4):331-348

Previous work on three‐dimensional shakedown analysis of cohesive‐frictional materials under moving surface loads has been entirely for isotropic materials. As a result, the effects of anisotropy, both elastic and plastic, of soil and pavement materials are ignored. This paper will, for the first time, develop three‐dimensional shakedown solutions to allow for the variation of elastic and plastic material properties with direction. Melan's lower‐bound shakedown theorem is used to derive shakedown solutions. In particular, a generalised, anisotropic Mohr–Coulomb yield criterion and cross‐anisotropic elastic stress fields are utilised to develop anisotropic shakedown solutions. It is found that shakedown solutions for anisotropic materials are dominated by Young's modulus ratio for the cases of subsurface failure and by shear modulus ratio for the cases of surface failure. Plastic anisotropy is mainly controlled by material cohesion ratio, the rise of which increases the shakedown limit until a maximum value is reached. The anisotropic shakedown limit varies with frictional coefficient, and the peak value may not occur for the case of normal loading only. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Zircon U–Pb ages and Hf isotope compositions of migmatites from the North Qinling terrane and their geological implications     

Q. Liu  Y.‐B. Wu  H. Wang  S. Gao  Z.‐W. Qin  X.‐C. Liu  S.‐H. Yang  H.‐J. Gong 《Journal of Metamorphic Geology》2014,32(2):177-193

Migmatites are predominant in the North Qinling (NQ) orogen, but their formation ages are poorly constrained. This paper presents a combined study of cathodoluminescence imaging, U–Pb age, trace element and Hf isotopes of zircon in migmatites from the NQ unit. In the migmatites, most zircon grains occur as new, homogeneous crystals, while some are present as overgrowth rims around inherited cores. Morphological and trace element features suggest that the zircon crystals are metamorphic and formed during partial melting. The inherited cores have oscillatory zoning and yield U–Pb ages of c. 900 Ma, representing their protolith ages. The early Neoproterozoic protoliths probably formed in an active continental margin, being a response to the assembly of the supercontinent Rodinia. The migmatite zircon yields Hf model ages of 1911 ± 20 to 990 ± 22 Ma, indicating that the protoliths were derived from reworking of Palaeoproterozoic to Neoproterozoic crustal materials. The anatexis zircon yields formation ages ranging from 455 ± 5 to 420 ± 4 Ma, with a peak at c. 435 Ma. Combined with previous results, we suggest that the migmatization of the NQ terrane occurred at c. 455–400 Ma. The migmatization was c. 50 Ma later than the c. 490 Ma ultra‐high‐P (UHP) metamorphism, indicating that they occurred in two independent tectonic events. By contrast, the migmatization was coeval with the granulite facies metamorphism and the granitic magmatism in the NQ unit, which collectively argue for their formation due to the northward subduction of the Shangdan Ocean. UHP rocks were distributed mainly along the northern margin and occasionally in the inner part of the NQ unit, indicating that they were exhumed along the northern edge and detached from the basement by the subsequent migmatization process.  相似文献   

Boundary frame contribution in coupled and uncoupled steel plate shear walls          下载免费PDF全文

Meng Wang  Daniel J. Borello  Larry A. Fahnestock 《地震工程与结构动力学》2017,46(14):2355-2380

The steel plate shear wall (SPSW) system is a robust option for earthquake resistance due to the strength, stiffness, ductility and energy dissipation that it provides. Although thin infill plates are efficient for resisting lateral loads, boundary frames that are proportioned based on capacity design requirements add significant structural weight that appears to be one of the factors limiting the use of the system in practice. An alternate configuration, the SPSW with coupling (SPSW‐WC), was explored recently as an option for increasing architectural flexibility while also improving overall system economy and seismic performance. The SPSW‐WC, which extensively employs flexural boundary frame contribution, has shown promise in analytical, numerical and experimental studies, but recent research on uncoupled SPSWs suggests that boundary frame contribution should not be considered for carrying seismic design shear. As a result, in the present study, boundary frame contribution in SPSWs was explored with detailed three‐dimensional finite element models, which were validated against large‐scale SPSW‐WC tests. Six‐story systems were considered, and the study matrix included single and double uncoupled SPSWs along with coupled SPSWs that had various degrees of coupling. Variations in design methodology were also explored. The modeling framework was employed to conduct static monotonic and cyclic pushover analyses and dynamic response history analysis. These analyses demonstrate the beneficial effect of coupling in SPSWs and illustrate the need to consider boundary frame contribution in design of coupled SPSWs. In addition, sharing design shear between the infill plate and the boundary frame is more generally shown to not be detrimental if this sharing is done in the design stage based on elastic analysis and the resulting boundary frame provides adequate secondary strength and stiffness following infill plate yielding. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

Variability in the vertical hyporheic water exchange affected by hydraulic conductivity and river morphology at a natural confluent meander bend          下载免费PDF全文

Jinxi Song  Guotao Zhang  Weize Wang  Qi Liu  Weiwei Jiang  Weiqiang Guo  Bin Tang  Haifeng Bai  Xinyi Dou 《水文研究》2017,31(19):3407-3420

River confluences and their associated tributaries are key morphodynamic nodes that play important roles in controlling hydraulic geometry and hyporheic water exchange in fluvial networks. However, the existing knowledge regarding hyporheic water exchange associated with river confluence morphology is relatively scarce. On January 14 and 15, 2016, the general hydraulic and morphological characteristics of the confluent meander bend (CMB) between the Juehe River and the Haohe River in the southern region of Xi'an City, Shaanxi Province, China, were investigated. The patterns and magnitudes of vertical hyporheic water exchange (VHWE) were estimated based on a one‐dimensional heat steady‐state model, whereas the sediment vertical hydraulic conductivity (K_v) was calculated via in situ permeameter tests. The results demonstrated that 6 hydrodynamic zones and their extensions were observed at the CMB during the test period. These zones were likely controlled by the obtuse junction angle and low momentum flux ratio, influencing the sediment grain size distribution of the CMB. The VHWE patterns at the test site during the test period mostly showed upwelling flow dominated by regional groundwater discharging into the river. The occurrence of longitudinal downwelling and upwelling patterns along the meander bend at the CMB was likely subjected to the comprehensive influences of the local sinuosity of the meander bend and regional groundwater discharge and finally formed regional and local flow paths. Additionally, in dominated upwelling areas, the change in VHWE magnitudes was nearly consistent with that in K_v values, and higher values of both variables generally occurred in erosional zones near the thalweg paths of the CMB, which were mostly made up of sand and gravel. This was potentially caused by the erosional and depositional processes subjected to confluence morphology. Furthermore, lower K_v values observed in downwelling areas at the CMB were attributed to sediment clogging caused by local downwelling flow. The confluence morphology and sediment K_v are thus likely the driving factors that cause local variations in the VHWE of fluvial systems.  相似文献   

Evaluation of Hydraulic Properties of Aquitards Using Earthquake‐Triggered Groundwater Variation          下载免费PDF全文

Shih‐Jung Wang  Kuo‐Chin Hsu  Chein‐Lee Wang  Wen‐Chi Lai  Liang‐Tzu Hsu 《Ground water》2017,55(5):747-756

The hydraulic properties of aquitards are not easily obtained because monitoring wells are usually installed in aquifers for groundwater resources management. Earthquake‐induced crust stress (strain) triggers groundwater level variations over a short period of time in a large area. These groundwater anomalies can be used to investigate aquifer systems. This study uses a poroelastic model to fit the postseismic variations of groundwater level triggered by the Chi‐Chi earthquake to evaluate the hydraulic properties of aquitards in the Jhoushuei River alluvial fan (JRAF), Taiwan. Six of the adopted eight wells with depths of 70 to 130 m showed good agreement with the recovery theory. The mean hydraulic conductivities (K) of the aquifers for the eight wells are 1.62 × 10^?4 to 9.06 × 10^?4 m/s, and the thicknesses are 18.8 to 46.1 m. The thicknesses of the aquitards are 11.3 to 42.0 m. Under the isotropic assumption for K, the estimated values of K for the aquitards are 3.0 × 10^?8 to 2.1 × 10^?6 m/s, corresponding to a silty medium. The results match the values obtained for the geological material of the drilling core and those reported in previous studies. The estimated values were combined with those given in previous studies to determine the distribution of K in the first two aquitards in the JRAF. The distribution patterns of the aquitards reflect the sedimentary environments and fit the geological material. The proposed technique can be used to evaluate the K value of aquitards using inverse methods. The inversion results can be used in hydrogeological analyses, contaminant modeling, and subsidence evaluation.  相似文献